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We propose a phase diagram for a cold polarized atomic Fermi gas with zero-range interaction. 
We identify four main phases in the plane of density and polarization: the superfluid phase, the 
normal phase, the gapless superfluid phase, and the modulated phase. We argue that there exists 
a Lifshitz point at the junction of the normal, the gapless superfluid, and the modulated phases, 
and a splitting point where the superfluid, the gapless superfluid, and the modulated phases meet. 
We show that the physics near the splitting point is universal and derive an effective field theory 
describing it. We also show that subregions with one and two Fermi surfaces exist within the normal 
and the gapless superfiuid phases. 

PACS numbers: 03.75.Ss 



I. INTRODUCTION 

The Fermi gas in the regime of large scattering length 
a [1-3] has attracted much interest due to its universal 
behavior. The regime can be achieved in atom traps by 
using the technique of Feshbach resonance [4-8]. Most 
attention is focused on systems consisting of two compo- 
nents of fermions (e.g., two spin components of a spin-i 
fermion) with equal number density. When the effective 
range tq is small compared to the interparticle distance 
n~^^^, where n is the total number density, many prop- 
erties of the system depend on n and a only through the 
dimensionless diluteness parameter 



When one varies n the system interpolates between 
the Bose-Einstein condensation (BEC) regime and the 
Bardeen-Cooper-Schrieffer (BCS) regime. For all values 
of K the ground state is believed to be a superfluid. 

In contrast, the case of unequal number density (or un- 
equal chemical potentials) of the two components is much 
less understood [9-11]. In the case of spin-i fermions 
one refers to a polarized gas. We follow this terminology, 
understanding "polarized" in the sense of asymmetry be- 
tween the two components. 

In this paper we propose a phase diagram for a po- 
larized Fermi gas in the whole range from the BEC to 
the BCS regime. Our proposal is summarized in Fig. 1. 
There are two special points on the phase diagram. Point 
S (the splitting point) is a point where phases 1, III and 
IV meet. Point L is a Lifshitz point where II, III and 
IV meet. The focus of this paper is on the point S. The 
physics in the vicinity of this point is of a long-distance, 
i.e., universal, nature and can be reliably studied within 
an effective field theory. 

Our proposal for the global structure of the phase di- 
agram is an educated guess anchored on the following 
reliable facts: the known phase structure in the BEC 
and BCS limits; the existence of the point S, which is at 
least a local minumum of the free energy; and the struc- 




FIG. 1: The proposed phase diagram in the plane of the di- 
luteness parameter k = —l/{nct') and the (scaled) polariza- 
tion chemical potential r) — II/Ah=o- I is the unpolarized 
BEC/BCS phase, II is the normal phase. III is the gapless 
superfluid phase, and IV is a region of Fulde-Ferrell-Larkin- 
Ovchinikov phases. The dashed line divides phases II and 
III into regions with one (on the left) and two (on the right) 
Fermi surfaces. Region IV must be divided into phases with 
different patterns of breaking of the rotational symmetry (not 
shown) . 

ture of the phase diagram around S, studied in detail in 
this paper. 

We devote the Appendix to the question of the global 
stability of the splitting point. We show that quantum 
Monte Carlo simulations [11] indicate that the superfluid 
state at point S is globally stable. We also argue that 
mean-field calculations of the phase diagram will likely 
miss the point S by significantly overestimating the size 
of the region occupied by phase II. 

II. QUALITATIVE DESCRIPTION OF THE 
PHASE DIAGRAM 

A. Axes of the phase diagram 

A particular system is characterized by three parame- 
ters: the scattering length a and the chemical potentials 
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of the two components /X| and Because of univer- 
sality (corresponding to rescaling invariancc a e^^a, 
jjLi e^^ iii) the whole phase diagram can be captured in 
a two-dimensional plot. We introduce the notation 

M = ^(MT+W), -H'=^(m-Mi)- (2) 

Then the parameter k on the horizontal axis is defined by 
Eq. (1) where n = n{ji, a) is the density of an unpolarized 
gas at chemical potential n and scattering length a. Thus 
K is the inverse diluteness parameter of an unpolarized 
system with chemical potential equal to the average of 
/X| and Hi and with the same scattering length a. 
The parameter r] on the vertical axis is defined as 



Here A(/i, a) is the gap in the fcrmion excitation spec- 
trum of the unpolarized gas. We divide if by A so that 
T] is dimensionless. This is also convenient because r] = 1 
is the line dividing phases I and III, as explained below. 

B. Qualitative understanding of the phase diagram 

The general structure of the phase diagram can be 
guessed by studying the BEC and BCS limits first. In the 
BEC limit k — oo, the gap in the fermionic excitation 
spectrum is close to the two-body binding energy. The 
excitations are fermionic quasiparticles carrying charge 
-|-1 and —1 with respect to the chemical potential H. 
When i7 < A or 77 < 1, the ground state is the same 
as for r] = since it is energetically unfavorable to cre- 
ate an unbound fermion on top of the BEC ground state 
(phase I). However, if 77 > 1, extra fermions will be cre- 
ated, and the system, for a range of 77, is a homogeneous 
mixture of bosonic bound states and fermions carrying 
one sign of spin (phase III). Finally, for sufficiently large r] 
all bosonic bound states disappear and the system is a 
completely polarized Fermi gas (phase II). 

The ground state in the BCS limit k — » +00 is also 
known. For r] < rji ~ l/-\/2, the ground state remains 
the BCS state with zero polarization (I). For rji < i] < 
0.754, the ground state is in one of the Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) states [12, 13] where Cooper 
pairs form with nonzero momentum and the superfluid 
order parameter varies in space (IV) . The precise spatial 
structure of the ground state is difficult to determine, 
but it is presumably crystalline. Finally for rj > 0.754 
the system goes to the Fermi liquid phase II. 

At small ?7, the unpolarized BEC phase can continue 
analytically into the unpolarized BCS phase. Indeed, in 
both regimes the U(l) symmetry associated with the con- 
servation of the total number of atoms is spontaneously 
broken. At large rj, the normal phase on the BEC side 
and that on the BCS side are essentially the same phase 
of a polarized Fermi gas. On both sides the U(l) sym- 
metry is restored. However, at intermediate rj the phases 




FIG. 2: Spectrum of fermionic quasiparticles at 77 = for 
K < Ko (left) and for k > ko (right). The minimum is located 
at p = for K < /to and p for k > Kq- 

in the two limits are qualitatively different: on the BCS 
side the translational invariance is spontaneously broken, 
while on the BEC side it is not, as already observed in 
Ref. [10]. 

Under the assumption that the phase diagram con- 
tains four phases which enter the region of intermedi- 
ate K and 77 — I from below, II from above. III from the 
left and IV from the right — there are two possibilities. 
The first is that there is a phase transition (coexistence) 
line separating the boson-fermion mixture phase III and 
the FFLO phase IV. The second is that III and IV do not 
coexist, but instead I and II do, which is what is found 
in mean-field theory. In the Appendix we argue that this 
is a likely artifact of the mean-field approximation. 

Within the first possibility, the minimal phase diagram 
therefore should look like Fig. 1. The line which separates 
the mixture phase III from the FFLO phase IV should 
start at a point S located on the boundary of the unpo- 
larized BEC/BCS phase I, and end at a point L on the 
boundary of the normal phase II. The point S will be 
called the splitting point, since, as we will show below, it 
is the point where the onset transition (from the unpo- 
larized phase I to the mixture phase III) splits into two 
first-order phase transitions. The point L, by the nature 
of the phases surrounding it, is a Lifshitz point [14]. 

In the next section we shall identify the splitting point 
S and study the phases surrounding it. We shall see that 
the physics around S matches very well with the global 
picture presented above. 



III. THE SPLITTING POINT 
A. Location of the splitting point 

Begin by considering the unpolarized gas, i.e., 77 = 0. 
The dispersion relation of a fermionic quasiparticle e{p) 
has a finite gap A in the whole range of k, from BEC 
to BCS. There is, however, a qualitative change in the 
location of the minimum of the dispersion curve as illus- 
trated in Fig. 2. In the BEC limit the dispersion curve is 
e{p) = I (2m) -\- A and achieves its minimum at p = 0. 
On the other hand, in the BCS limit the dispersion curve 
has a minimum at p 7^ (near the Fermi momentum). 

One concludes that, as the system crosses over from 
the BEC to the BCS regime, at a certain point n = kq 
the minimum of the dispersion curve shifts from p = 
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to p 7^ 0. In the mean-field theory this occurs when 
the chemical potential vanishes [2] ; we shall only assume 
that this transition occurs smoothly. Around this n the 
dispersion curve has the following behavior near p = 0: 



= eo - ap^ + Pp^ 



(4) 



where a changes sign at k = kq: a = q:o(k— kq), do > 0. 
Note that A = eo - 0{a) aV(4/3)- 

At = 0, the change in the location of the minimum 
of the dispersion curve by itself does not signal any phase 
transition. However, as we shall see shortly, the splitting 
point S is located at k = kq and rj ~ 1. 

The value of Kq cannot be found from effective theory. 
One can determine the sign of n from the shape of the 
dispersion curve at infinite scattering length [n = 0). 
According to a recent quantum Monte-Carlo calcula- 
tion [11], at K = the dispersion is BCS-like, i.e., has 
a minimum at a nonzero value of p. From this informa- 
tion we conclude that kq < 0, i.e., on the left (BEC) side 
of the phase diagram in Fig. 1. 



B. Effective field theory near the splitting point 

Let us take the system at some n close to Kq and in- 
crease r]. When r] is close to 1, there are fermion quasipar- 
ticles whose energy e = e(p) — H is small. The momenta 
of these quasiparticles are also small, around the mini- 
mum of the dispersion curves. At small number density, 
these quasiparticles are weakly interacting. The whole 
problem is therefore treatable by using a low-energy ef- 
fective field theory of quasiparticles (supcrfluid phonons 
and the extra fermions) despite the fact that the original 
(undressed) particles are strongly coupled. 

The effective Lagrangian is constrained by two U(l) 
symmetries. One of them corresponds to the conser- 
vation of the total number of atoms. The supcrfluid 
mode ip, which we normalize to be half of the phase 
of the Cooper pair, transforms under this symmetry as 
(fi ^ (fi + a. The fermion field ip can always be chosen 
to be neutral under this symmetry by multiplying it by 
an appropriate power of e*'^. Another U(l) symmetry 
corresponds to the conservation of the difference of the 
numbers of the two components, or polarization. Un- 
der this symmetry ip — > ipe'^^ , while (f is invariant. To 
lowest nontrivial order in (p and ijj and derivatives the 
Lagrangian obeying these symmetries is given by 



£ = 



The low-energy parameters ft, f and g are not con- 
strained by the U(l) symmetries. Their values are not 
essential for the discussion of the phase diagram below. 
However, for the purpose of further applications, we point 
out that they are constrained by the Galilean invariance. 
Using the results of Refs. [15, 16], one can show that 
ff = dn/dfi, = n/m, and 5 = 1. 



Note that the interaction term ^jj'^ijj \lp + {dipY / {2m)\ 
is allowed, but is negligible when the number density of 
fermion quasiparticles is small, ■^^■^ <C n. The fermion 

self-interaction {ij)^ d t/j)"^ is also irrelevant. 

C. Phases near the splitting point 

When a > the ground state may carry nonzero spa- 
tial gradient of ip or, in other words, nonzero supcrfluid 
velocity Vg = dip/m. This can be seen, e.g., from the fact 
that a state with an arbitrarily small density of fermions 
and Vs = has a negative superfluid density due to the 
divergent density of states at the Fermi surface [17]. Our 
task is to find, at each value of a (close to 0), the mini- 
mum of the free energy as a function of Vg . Using Eqs. (4) 
and (5) we find the fermion dispersion relation in the 
presence of a uniform superfluid flow with velocity Vg : 



ev{p) = eo - ap^ +f3p'^ 



+ Vs P- H. 



(6) 



All levels with ey{p) < are filled. The total free cnc;rgy 
receives contributions from the superfluid flow and from 
the filled fermionic energy levels. 



d^p 

(27r)3 



ey{p)0{-e,{p)). (7) 



Here p = mn, and 9{x) is the step function. Performing 
the integration, we find the free energy. 



4(157r2)4 p3p' 



(8) 



where we rescaled Vg and {H — A) by introducing dimen- 
sionless variables 



h = 2(157r^)^^(ff- A), 



a" 



(9) 
(10) 



The shape of F{vg) at a given H is captured by the func- 
tion fh{x), 



fh{x) 



X 



(h + xf'^6{h + x)-{h- xf'^6{h - x) 

(11) 

at the corresponding h, which is shown in Fig. 3. For 
h < h\ !v —0.067 the absolute minimum of / is located 

at a; = 0. When hi < h < /),2 « 0.502 the minimum 
switches to a; 7^ 0. When /i > /12, the minimum switches 
back to X = 0. 

The change in the location of the minimum of / trans- 
lates into two first-order phase transitions occurring at 



r?i,2 = 1 + 



hi:. 



2(157r2)2 p2^4^ • 



(12) 



The following picture emerges. At fixed k > kq, as we 

increase rj from the ground state is an impolarized su- 
perfluid (phase I) for rj < rji, until a phase transition at 



Ux) - /fc(0) h = -0.078 




V = V3 



FIG. 3: Function fh{x) at three representative values of h. 



77 = 771. Note that, since hi < 0, the phase transition 
occurs while H is still below the gap: Hi < A. 

For Tji < rj < ri2 the ground state has nonzero super- 
fluid velocity (phase IV). It also contains fermionic 
quasiparticles filling an asymmetric region in the phase 
space [Figs. 4 and 5(a)]. The currents carried by the su- 
perfluid flow and by fermion quasiparticles cancel each 
other so that all currents (both total number and polar- 
ization) in the ground state are zero. 














' Px 



FIG. 4: Occupied fermion quasiparticle modes in phase IV 
near point S occupy a thin spherical cap with curvature radius 
p = \/a/(2/3). The angle 9 increases from 80.4° at 77 = r;i + 
to 127.3° at 77 = 772 - 0. 

When rj > rj2 the superfluid velocity switches back to 
zero. However, now > A, so there is a flnite density 
of fermionic quasiparticles (phase III). At 77 = 7/2 -I- 
the fermions flU a thin spherical shell centered at p = 
[Fig. 5(b)]. As 77 increases, the shell thickens and at 



fo 
A 



1 



a 

4/3A 



(13) 



or H = eo, the shell turns into a solid ball [Fig. 5(c)]. 
This occurs at the dashed line on Fig. 1. 



IV. GLOBAL TOPOLOGY OF THE PHASE 
DIAGRAM AND CONCLUDING REMARKS 

In the region 771 < 7y < 7^2 where Vs 7^ 0, the condensate 
function varies as e**^ with q = 2mVs 0. This has the 




FIG. 5: Occupied quasifermion modes {p^ = cross section) 
near point S in phase IV (a) and in phase III (b),(c). 



same structure as the Fulde-Ferrell phase [12] in the BCS 
regime. Therefore, it is natural to expect that the region 
ryi < 77 < 772 near S is continuously connected to the 
FFLO region on the BCS side, as in Fig. 1. 

In the BCS limit, the ground state is expected to be a 
crystal where A(a;) is a superposition of several plane 
waves [18]. However, A{x) has to be a single plane 
wave near S. Indeed, Vs is small (vs — > as one ap- 
proaces S), while the healing length is flnite at S. If one 
tried to superimpose two plane waves to make a stand- 
ing wave cos{qx), the absolute value of the condensate 
would vary on length scales much larger than the healing 
length, costing a large amount of energy. Therefore, one 
concludes that the FFLO region of the phase diagram 
is further subdivided into smaller regions with different 
patterns of breaking of rotational and translational sym- 
metries. 

The line separating the gapless superfluid phase HI and 
the phase with nonzero Vs IV should end in a Lifshitz 
point L. It is known that in the vicinity of a Lifshitz 
point, the transition between the phase with spatially 
varying order parameter and that with constant order 
parameter can be second order when the broken sym- 
metry is continuous [14]. If this is the case, then the 
flrst-order line rj = 7/2 (k) should become second order at 
some value of k before the point L. Another alternative 
is that the entire line from S" to L is flrst order. 

In the region between 77 = 772 (k) and rj — 773 (k), but be- 
low the transition to the normal state (part of region III 
below the dashed curve on Fig. 1), the spectrum of low- 
energy excitations consists of the superfluid boson and 
the fermionic quasiparticles with momenta around two 
Fermi surfaces [the inner and the outer surfaces of the 
spherical shell in Fig. 5(b)]. This is the Sarma phase [19], 
which is unstable in the BCS limit but, as we found, is 
stable near the splitting point. This is interesting since 
previous examples of stabilized Sarma, or breached pair- 
ing, phase require quite special conditions [20] . 

Figure 1 describes the system at zero temperature T. 
At small r 7^ the phase diagram undergoes obvious 
modifications. Since the phases I and III have the same 
symmetry, the transition between them becomes an an- 
alytic crossover at T 7^ 0. Similarly, the dashed line on 
Fig. 1 should turn into a crossover at T ^ 0. On the 
contrary, the phases II and III must remain separated by 
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a phase transition because the U(l) symmetry is sponta- 
neously broken in III and restored in II. Likewise, since 
the rotational symmetry undergoes spontaneous break- 
down in phase IV, this phase must remain separated from 
all other phases by a phase transition. 

The results obtained in this paper should be interest- 
ing for other physical systems with asymmetric fermion 
pairing, for example for the cores of neutron stars [21] 
or QCD at finite baryon and isospin chemical potentials 
[22] . We defer this investigation to future work. 
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APPENDIX: GLOBAL STABILITY OF THE 
SPLITTING POINT 

Near the splitting point we determined the minima 
of the free energy with respect to small perturbations 
around the superfluid state. By its nature, the effective 
field theory wc applied cannot assert that the minima we 
study are global, i.e.. stable with respect to large devi- 
ations from the state we considered. We must therefore 
employ different methods to assess the global stability 
of the point S. It is an important question since, as we 
show below, a mean-field approximation will likely show 
the normal, not superfluid, phase to be the globally stable 
state at the point S. Since the mean-field approximation 
is not a controllable approximation in the relevant regime 
|k| < 1, other approaches are necessary to determine the 
global structure of the phase diagram. 

Since the point S lies on the line rj = H/A = 1, the 
relevant comparision is between the pressure of the super- 
fluid state Psf(M) H) at H = A and the pressure of the 
normal state P^ip, H) at the same H. We know that the 
normal state wins in the BCS limit, and the superfluid 
state wins in the BEG limit. Thus, there is a crossover 
at some value of Kc 



PsfIMj A) = Pn(/x, A) at K = Kc 



(A.l) 



Assuming there is only one such crossover point, the sta- 
bility of the point S requires that kq < Kc, i-e., that the 
splitting point S occurs while the superfluid state is still 
globally stable, as in Fig. 1. 

Although no direct determination of Kc has been made 
to date, Carlson and Rcddy [11] have performed the pres- 
sure comparision at the unitarity point k = 0, using a 
quantum Monte Garlo method. We shall now show that 
these results imply k^. w 0, as well as that kq < 0, and 
thus that the splitting point is globally stable. 



At the point k = the only dimensionful parameters 
are the chemical potentials /i and H. The pressure of 
the superfluid phase is the same for all H e (0, A) due 
to the gap, and, by dimensionality, is proportional to the 
pressure of a free two-component unpolarized Fermi gas 
at the same chemical potential. 



Psf(/x,A) = Psf(/x,0) 



Pfree(M) 
^3/2 



(A.2) 



where ^ is a universal dimcnsionlcss constant, conven- 
tionally defined as the ratio of the energy of the gas at 
unitarity and the energy of the noninteracting gas at the 
same density. Pfrcc(Ai) = 2{2m)^/^ / {Ibn'^) fi^/^ . 

To find the pressure of the normal phase, wc first no- 
tice that Hi = fj, — A < since A > /i at unitarity 
(both in mean-field and in quantum Monte Carlo meth- 
ods). Therefore the normal phase is a completely polar- 
ized (single-component — hence a factor 1/2 below) Fermi 
gas with Hi = h + A, whose pressure is 



Pn(m,A) = iPfree(MT) = ^ ' ^ 



5/2 



Pfree(M) (A.3) 



The ratio of the pressures of the superfluid and the 
normal phases is therefore expressed in terms of two di- 
mensionless quantities, ^ and A/fi, 



Psf(m,A) 
Pn(m,A) 



2 



1 + ^ 



-5/2 



at K = 0. 



(A.4) 



The Monte Garlo simulation of Ref. [11], finds ^ « 0.42 
and A/h ^ 1-2. Substituting into Eq. (A.4) one finds the 
ratio to be consistent with 1 within errors [11]. Compar- 
ing to Eq. (A.l) we conclude that Kc ~ 0. 

Carlson and Rcddy [11] also found that at k = the 
fermion dispersion curve has a miminum at nonzero mo- 
mentum, which means that the point k = is on the BCS 
side of the splitting point, i.e., kq < 0. Since kq < « Kc, 
we can conclude that the superfluid phase is still globally 
stable at the splitting point. 

On the other hand, if we substituted into Eq. (A.4) 
the mean-field values ^ ~ 0.59, A/fi « 1.16 [23], we 
would find the ratio to be 0.64, i.e., the supcrfiiud phase 
would appear to lose by a large margin at k = 0. There- 
fore it is likely that in a mean-field calculation of the 
phase diagram the splitting point will be overshadowed 
by the normal phase. Instead of Fig. 1, in the regime 
|«;| < 1 one would then observe a single first-order transi- 
tion (phase co-existence) line separating phases I and II. 
Qualitatively, it should be expected that a variational cal- 
culation overestimates the energy of the superfluid phase 
(i.e., the value of ^). It is an interesting problem to see 
if the mean- field theory [9, 10] can be improved in the 
regime |k| < 1 to describe the phase diagram correctly. 

A calculation of the phase diagram using a two-channel 
Hamiltonian was reported recently [24]. It relies on a 
mean-field approximation which becomes controllable in 
the weak channel coupling (narrow resonance) limit. It 
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should be borne in mind that in this hniit the efFectivc 
range of the interaction diverges: Tq ^ This 



hmit [25] is therefore opposite to the zero-range Umit we 
consider in this paper. 
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